The fate of (phantom) dark energy universe with string curvature corrections 
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We study the evolution of (phantom) dark energy universe by taking into account the higher- 
order string corrections to Einstein-Hilbert action with fixed dilaton and modulus fields. While the 
presence of a cosmological constant gives stable de-Sitter fixed points in the cases of heterotic and 
bosonic strings, no stable de-Sitter solutions exist when a phantom fluid is present. We find that 
the universe can exhibit a Big Crunch singularity with a finite time for type II string, whereas it 
reaches a Big Rip singularity for heterotic and bosonic strings. Thus the fate of dark energy universe 
crucially depends upon the type of string theory under consideration. 
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I. INTRODUCTION 

Recent observations suggest that the current universe 
is dominated by dark energy responsible for an acceler- 
ated expansion The equation of state parameter w 
for dark energy lies in a narrow region around w = — 1 
and may even be smaller than —1 Q. When w is less 
than —1, dubbed as phantom dark energy, the universe 
ends up with a Big Rip singularity 0, which is char- 
acterized by the divergence of curvature of the universe 
after a finite interval of time (see Refs. jHEJ)- 

The energy scale may grow up to the Planck scale in 
the presence of phantom dark energy. This means that 
higher-order curvature or quantum corrections can be im- 
portant around the Big Rip. For example, quantum cor- 
rections coming from conformal anomaly are taken into 
account in Refs. |jj for dark energy dynamics. It was 
found that such corrections can moderate the singularity 
by providing a negative energy density 8] . Thus it is im- 
portant to implement quantum effects in order to predict 
the final fate of the universe. 

In low-energy effective string theory there exist higher- 
curvature corrections to the usual scalar curvature term. 
The leading quadratic correction corresponds to the 
product of dilaton/modulus and Gauss-Bonnet (GB) cur- 
vature invariant ||| ■ The GB term is topologically invari- 
ant in four dimensions and hence does not contribute to 
dynamical equations of motion if the dilaton/modulus 
field is constant ^(| . Meanwhile it affects the cosmologi- 
cal dynamics in presence of dynamically evolving dilaton 
and modulus fields. The possible effects of the GB term 
for early universe cosm ology an d black hole physics were 
investigated in Refs. [llL Il2| . Lately the GB correc- 
tion was applied to the study of cosmological dynamics 
of dark energy . 

When the dilaton and mudulus are fixed, it is impor- 
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tant to implement third and next-order string curvature 
corrections 0] ■ This can change the resulting cosmolog- 
ical dynamics drastically as it happens in the context of 
inflation 0] and black holes 0| . The goal of the present 
paper is to study the effect of next-to-leading order string 
corrections to the cosmological dynamics around the Big 
Rip singularity with an assumption that the dilaton and 
the modulus are stabilized. We would also investigate 
the existence and the stability of de-Sitter solutions in 
the presence of a cosmological constant. We shall con- 
sider three types of string corrections and study the fate 
of the universe accordingly. 



II. EVOLUTION EQUATIONS 

Let us consider the Einstein-Hilbert action in low- 
energy effective string theory: 



S= / d D Xy/=jj[R+C c + ...] , (1) 



where R is the scalar curvature and C c is the string cor- 
rection which is given by |lfjj | 

C c = c lC /e-^/; 2 + c 2 a' 2 e~ 4 *C 3 + c 3 a ,3 e^C 4 , (2) 

where a' is the string expansion parameter, (j) is the dila- 
ton field, and 

C 2 = Q 2 , (3) 

C 3 = 2n 3 +R^Rf p R^, (4) 

£4 = £41 — 5 H C i2 — 5 (5) 
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Here one has 8h(b) — 1 f° r heterotic (bosonic) string and 
zero otherwise. The Gauss-bonnet term, J7 2 , as well as 
the Euler density, O3, does not contribute to the back- 
ground equation of motion for D = 4 unless the dila- 
ton is dynamically evolving. The coefficien ts ( cj , c 2 , C3) 
are different depending on string theories |10|. We have 
(ci,c 2 ,c 3 ) = (0,0, 1/8), (1/8, 0,1/8), (1/4, 1/48, 1/8) for 
type II, heterotic, and bosonic strings, respectively. In 
the case of type II string with D = 4, for example, only 
the £41 term affects the dynamical evolution of the sys- 
tem. 

We shall consider the flat Friedmann-Robertson- 
Walker metric with a lapse function N(i): 



ds 2 = -N{t) 2 dt 2 + a{tf ^T(do 



,i\2 



(11) 



where d = D — 1. The Ricci tensors under this metric 
are given in the Appendix. In what follows we shall con- 
sider the case of D = 4 under the assumption that the 
modulus field which corresponds to the radius of extra di- 
mensions is stabilized after the compactification to four 
dimensions. Then we find 



£ 3 - J^(H +I)-—HI 
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where H = a /a is the Hubble rate and / is defined 
by I = H 2 + H . It should be noted that, in the case 
of de-Sitter space time, the expressions ifT2 |l -lfT5 jl reduce 
to their counterparts given by Eqs. (|27|I - H30|I in the Ap- 
pendix. 

We shall implement the contribution of a barotropic 
perfect fluid to the action (JJJ. The equation of state 
parameter, w m = p m / Pm, is assumed to be constant. 
Our main interest is to study the final fate of universe 
filled with a phantom- type fluid (w m < — 1). In this 
case the universe eventually reaches a Big Rip singular- 
ity 0] with a divergent Hubble rate in the absence of 
higher-curvature terms. We are interested in the effect 
of string curvature corrections to the cosmological evolu- 
tion around the Big Rip. 

Varying the action (JIJ with respect to N, we find 



6H 2 = Pc 



where 



dt\dN J 



ON dN 



(16) 
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p m is the energy density of the barotropic fluid, satisfying 



Pn 
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u m JHrn 



= 0. 



(18) 



In what follows we shall consider the case with a fixed 
dilaton. Then we have two dynamical equations i|16|) 
and (|18fl for our system. We note that the variation of 
the action (JIJ in terms of the scale factor a gives rise to 
another equation, but this can be derived from Eqs. Ijlfijl 
and (|18f) by taking a derivative with respect to t. From 
Eq. I|17l) we find that the energy density p c for type II & 
heterotic strings is 

p c = B[a$H 8 + a c I 4 + a 4 H 4 I 2 + a 2 H 2 I 3 + a 6 H 6 I 

-J{a 5 H 5 + ai HI 2 + a 3 H 3 I)}, (19) 

where J = H + 3HH+H 3 . One has B = 6c 3 a' 3 e- 6 ^C(3), 
a§ = — 21, a c = — 3, a4 = —12, 02 = 4, ciq = —24, 05 = 
— 8, ai = — 12, 0,3 — —24 for type II string, and B = 
6c 3 a' 3 e- 6 <t>, a 8 = -21((3) + 35, a c = -3((3) + 15, a 4 = 
-12C(3) - 6,a 2 = 4C(3) + 20, a 6 = -24C(3) + 12, a 5 = 
-8C(3) + 4, 01 = -12C(3)+60, a 3 = -24((3) for heterotic 
string. In the case of bosonic strings we have 

p c = A(5H e + 2I 3 -QHIJ)+B[(-2l((3)+210)H s 
+ (-3C(3) + 90)/ 4 - (12C(3) + A8)H 4 I 2 
+(4C(3) + \2Q)H 2 I 3 + (-24C(3) + 96)i? 6 / 
+J{(8((3) - 32)i? 5 + (12C(3) - 360)HI 2 
+24C(3) J ff 3 /}] . (20) 

where A = 24c 2 a' 2 and B = Qc 3 a' 3 e~ 64 '. 



III. THE FATE OF DARK ENERGY UNIVERSE 

In this section we study the cosmological evolution in 
dark energy universe for the above three classes of string 
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curvature corrections. Our main interest is the universe 
dominated by phantom dark energy, but we consider the 
case of cosmological constant as well. We first study the 
effects of type II and heterotic corrections and then pro- 
ceed to the bosonic correction. 



A. Type II and heterotic strings 

For the analysis of dynamics to follow, it would be 
convenient for us to cast Eqs. (|16[) and (|18|l with the 
correction term (I19|) in the form: 

x = y, (21) 
y = [B{a 8 x 8 + a c I 4 + a A x A l 2 + a 2 x 2 I 3 + a 6 x 6 I 

-(3xy + x 3 )Q+z-6x 2 ]/(BZ), (22) 
z = -3(1 + w m )xz , (23) 

where x = H, y = H, z = p m , I = x 2 + y and £ = 
05a; 5 + aixl 2 + a^x 3 I. Here we shall consider the case 
with w m 7^ —1. By setting x = 0, y = and z — 0, we 
find the following de-Sitter fixed point: 

/ 6 \ 1/6 

x c = , Vc = 0, 2 C = 0, (24) 

where D = a c + a 2 + a A + ae + a 8 — &i — «3 — 15- Since 
D < for both type II and heterotic strings, we do not 
have de-Sitter fixed points. 

When a cosmological constant A is present instead of 
p m , corresponding to the equation of state w rn = —1, we 
find from Eq. I|16fl that there exists one de-Sitter solution 
which satisfies A = 6 if 2 — BDH 8 . One can study the sta- 
bility of this solution by considering small perturbations 
6x and Sy about the fixed point. We evaluate two eigen- 
values for the matrix of perturbations using the method 
in Ref. For the type II correction we find that one 
eigenvalue is positive while another is negative, thereby 
indicating that the de-Sitter solution is not stable. Mean- 
while in the heterotic case the de-Sitter solution is either 
a stable spiral (for smaller values of A, see Fig. ^) or a 
stable node (for larger values of A). We note, however, 
that this stable solution disappears for the equation of 
state with w m ^= — 1. 

In order to understand the fate of the universe which 
is dominated by a phantom-type fluid, we solve au- 
tonomous equations (|2~TJ> - lj2"3l numerically. Figure El 
shows the evolution of the Hubble rate for a phantom- 
type fluid with w m — —1.5 in the presence of string cur- 
vature corrections. We find that in the type II case the 
Hubble rate begins to decrease because of the presence of 
string corrections and it eventually diverges toward — 00 
after crossing H — 0. Thus the fate of the universe is 
characterized by a Big Crunch rather than a Big Rip. 
Meanwhile in the heterotic case the Hubble rate contin- 
ues to increase and diverges after a finite interval of time, 
see Fig. Thus the Big Rip singularity is inevitable even 
when the heterotic string correction is present. 
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FIG. 1: The phase portrait for heterotic string in the case 
of p m = A = 1 for several different initial conditions. The 
stable fixed point x c = H c = 0.408 corresponds to a de-Sitter 
solution which is a stable spiral. 
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FIG. 2: The evolution of the Hubble rate H in the presence 
of string curvature corrections and a phantom fluid with an 
equation of state: w m = —1.5. For the type II correction, the 
solution approaches H = — 00 by crossing H — 0, whereas in 
the heterotic and bosonic cases the Hubble rate grows toward 
infinity. 

We also studied cosmological evolution for several dif- 
ferent values of w m and for different initial conditions 
of the Hubble rate. For the type II correction we find 
that the solutions approach the Big Crunch singularity 
for w m < —1-2, whereas they tend to approach the Big 
Rip singularity for —1.2 < w m < —1. Thus the fate of 
the universe depends upon the equation of state for phan- 
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torn dark energy. For the heterotic case we find that the 
solutions reach to the Big Rip singularity independent of 
the values of w m (< — 1) and initial conditions of H. 

B. Bosonic string 

As demonstrated above, the type II and the heterotic 
string models do not exhibit de-Sitter solutions for w m ^ 
— 1. However, in the case of bosonic string, there exists 
a de-Sitter fixed point which satisfies the relation 

B [76 - 12C(3)] H 6 + AH 4 -6 = 0. (25) 

For example one has H = 0.709 for A = 1/2 and 
B = 3/4. By considering small perturbations around this 
solution and evaluating three eigenvalues of the 3x3 ma- 
trix for the system given by Eqs. (|21|l - 123|) , we find that 
two of the eigenvalues are positive for w m < — 1 and one 
of them is positive for w rn > — 1. Therefore the de-Sitter 
solution characterized by Eq. I|25|l does not correspond 
to a stable attractor for w m ^ — 1. 

In Fig.[2]we plot the evolution of the Hubble rate with 
bosonic string corrections in the presence of a phantom- 
type fluid with w m = —1.5. The Hubble rate continues 
to grow and diverges with a finite time as in the case of 
heterotic string. We also run our numerical code for dif- 
ferent values of w m and find that the solutions approach 
the Big Rip singularity for w m < — 1. 

When a cosmological constant A is present instead of 
p m , de-Sitter solutions satisfy 

A = 6H 2 - AH 6 - B [76 - 12C(3)] H s . (26) 

There exist two solutions for this equation provided that 
A ranges in the region < A < /(Hm) = — 
AH Q M - B [76 - 12C(3)] J ff| f , where H M is the solution 
for ZAHij + 4B[76 - 12C(3)]Jl£ f = 6. For example 
H M = 0.56 and f{H M ) = 1.42 for A = 1/2 and B = 3/4. 
When A = 1, one has two de-Sitter solutions character- 
ized by H = 0.417 and H = 0.652. We evaluate two 
eigenvalues of the 2x2 matrix for perturbations Sx and 
Sy around the fixed points. We find a complex conjugate 
pair of eigenvalues with negative real part for H = 0.417 
making the fixed point a stable spiral. As for the second 
second critical point corresponding to H = 0.652, one 
of the eigenvalues turns out to be positive, which means 
that the de-Sitter solution is unstable— a saddle in this 
case (see Fig. 01. 

IV. SUMMARY 

In this paper we studied the effect of higher-curvature 
corrections in low-energy effective string theory on the 
cosmological dynamics in the presence of dark energy 
fluid. Since the existence of a phantom fluid leads to the 
growth of the Hubble rate, the energy scale of the uni- 
verse may reach the Planck scale in future. This means 
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FIG. 3: The phase portrait for bosonic string in the presence 
of a cosmological constant (p m sA = l) with A = 1/2 and 
B — 3/4. There exist two de-Sitter fixed points. The point A 
[x c = 0.417] corresponds to a stable spiral, whereas the point 
B [x c = 0.652] is a saddle. 



that string curvature corrections can be very important 
to determine the dynamical evolution of the universe. 

We have considered string corrections up to quartic in 
curvatures for three type of string theories- (i) type II, 
(ii) heterotic, and (iii) bosonic strings. In our analysis the 
contribution of the Gauss-Bonnet term does not affect 
the cosmological dynamics in D — 4 dimensions, since 
the dilaton and the modulus fields are fixed. For the fluid 
with an equation of state characterized by w m ^ —1, we 
find that de-Sitter solutions do not exist for type II and 
heterotic string corrections. There is a de-Sitter solution 
for the bosonic string even for w m 7^—1, but this is found 
to be unstable. When the equation of state for fluid is 
that of a cosmological constant (w m = —1), we find that 
stable de-Sitter solutions exist for heterotic and bosonic 
strings. 

We ran our numerical code to study the effect of string 
corrections around the Big Rip when a phantom fluid is 
present. In the type II case we found that the solutions 
approach the Big Crunch singularity (H — > — oo) after 
crossing H = when the equation of state for dark en- 
ergy is w m < —1.2. Meanwhile the Hubble rate diverges 
toward H — > +oo with a finite time for heterotic and 
bosonic corrections, which implies that the Big Rip sin- 
gularity is difficult to be avoided in these cases. The 
divergent behavior of the Hubble rate for w m < — 1 is 
associated with the fact that there are neither stable de- 
Sitter nor stable Minkowski attractors for the types of 
the corrections we considered. 

In the present letter we restricted our attention 
to purely geometrical effects assuming that non- 



perturbative potentials may arise allowing to freeze dila- 
ton and modulus fields. Compactifications from higher 
dimensions to 4-dimensional space time result in residual 
modulus fields which are related to the radii of internal 
space. In general a modulus field is dynamical and inter- 
acts with higher-order curvature terms. The same is true 
for the dilaton field which is related to string coupling g s . 
These may give rise to non-trivial effects, for instance, 
even the GB curvature invariant which is purely topolog- 
ical in 4 dimensions, does not vanish in the presence of 
dynamical modulus and dilaton fields. Such a scenario 
would have important implications for future evolution 
of the dark energy universe. Very recently, cosmologi- 
cal dynamics based upon effective string theory action 
was investigated with dynamically evolving modulus and 
dilaton fields in the presence of second-order curvature 
corrections ^3 ( see R e f- E3 on the related theme). It 
was demonstrated that the second-order curvature cor- 
rection to Einstcin-Hilbert action can significantly mod- 
ify the structure of future singularities in dark energy 
universe. It is therefore important, though technically 
cumbersome, to extend the analysis of the present letter 
to the case of dynamical modulus/dilaton fields. 

A comment is also in order about the distinct features 
that different string models exhibit in cosmological dy- 
namics. With fixed modulus/dilaton, the evolution of 
phantom dark energy universe is clearly distinguished de- 
pending upon string models, namely, the fate of such a 
universe is Big Crunch for the type II superstring whereas 
it is a Big Rip for bosonic and heterotic strings. This 
distinction mainly comes from the difference of the coef- 
ficients in Eqs. ljT§|) and (J2SJ- In the type II case p c be- 
comes negative, which counteracts the energy density of 
the phantom fluid. This property is crucially important 
to avoid the Big Rip singularity as pointed out in Ref. Q . 
We note that this behavior also appears in the presence 
of a dynamical modulus field with second-order string 
corrections . It is really of interest to investigate how 
the final fate of the universe is changed when dynami- 
cal modulus/dilaton fields couples to third/fourth-order 
string curvature corrections. We hope to address this 
issue in future work. 

In addition, the higher-order curvature contributions 
used in our description have inbuilt ambiguities related 
to particular metric redefinitions. It would be important 
to investigate whether or not these ambiguities can lead 
to different fate of cosmological evolution. 



where fi — 1,2, ... d. By using the formula 
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we find that the non-zero components of Riemann tensors 
are 
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Noting the relation Rmn — Rmln> we obtain 
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We also find 
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These relations are used to evaluate the correction term 
C c . 

We should note that the corrections can be eas- 
ily computed in case of de-Sitter symmetry [R a bcd = 
A (g ac gbd - 9ad9bc)], as 
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APPENDIX: Calculation of curvature tensors 



C 3 = 4A 3 d(d+l), (27) 
£ 4 i = -3C(3)A 4 (d- l)d(d + I), (28) 

A2 = ^A i (d-7)d 2 (d+l), (29) 

£ 43 = -4A 4 d(d + 1)[1- d(d- 9)] . (30) 

These coincide with Eqs. l(T2 jl -(fT5 }l by setting N — 1 and 
H 2 = A. 



For the metric (|II|I the non-zero components of 
Christoffel symbols are 
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